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Abstract 

In this paper we review studies on spin-dependent transport in systems contain- 
ing ferromagnetic nanoparticles. In a tunnel junction with a nanometer-scale- 
island, the charging effect leads to an electric current blockade phenomenon in 
which a single electron charge plays a significant role in electron transport, re- 
sulting in single-electron tunneling (SET) properties such as Coulomb blockade 
and Coulomb staircase. In a tunnel junction with a ferromagnetic nano-island 
and electrode, it was expected that the interplay of spin-dependent tunneling 
(SDT) and single-electron tunneling (SET), i.e., spin-dependent single-electron 
tunneling (SD-SET), would give rise to remarkable tunnel magnetoresistance 
(TMR) phenomena. We investigated magnetotransport properties in both se- 
quential tunneling and cotunneling regimes of SET and found the enhancement 
and oscillation of TMR. The self-assembled ferromagnetic nanoparticles we have 
employed in this study consisted of a Co-Al-0 granular film with cobalt nanopar- 
ticles embedded in an Al-0 insulating matrix. A C036AI22O42 film prepared by 
a reactive sputtering method produced a TMR ratio reaching 10 % and su- 
perparamagnetic behavior at room temperature. The TMR ratio exhibited an 
anomalous increase at low temperatures but no indication of change with bias 
voltage. In Ch. 4, we show that the anomalous increase of the MR provided ev- 
idence for higher-order tunneling (cotunneling) between large granules through 
intervening small granules. We emphasize that the existence of higher-order 
tunneling is a natural consequence of the granular structure, since broad distri- 
bution of granule size is an intrinsic property of granular systems. In Ch. 5, we 
concentrate on SD-SET properties in sequential tunneling regimes. We fabri- 
cated two types of device structures with Co- Al-0 film using focused ion-beam 
milling or electron-beam lithography techniques. One had a granular nanobridge 
structure: point-shaped electrodes separated by a very narrow lateral gap filled 
with the Co- Al-0 granular film. The other had a current-perpendicular-to-plane 
(CPP) geometry structure: a thin Co-Al-0 granular film sandwiched by ferro- 
magnetic electrodes with the current flowing in the direction perpendicular to 
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the film plane through a few Co particles. We found the enhancement and os- 
cillation of TMR due to spin-dependent SET in sequential tunneling regimes. 
In Ch. 6, we report experimental evidence of a spin accumulation effect in Co 
nanoparticles leading to the oscillation of TMR with alternate sign changes. 
Furthermore, we discovered that the spin relaxation time in the nanoparticles 
is unprecedentedly enhanced up to the order of more than hundreds of nanosec- 
onds, compared to that evaluated from the spin-diffusion length of ferromagnetic 
layers in previous CPP-GMR studies, i.e., the order of tens of picoseconds. 
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1. Introduction 

Ever Since the giant magnetoresistance (GMR) effect was discovered in 
1988 [1], magnctotransport phenomena in magnetic multilayers with precise layer- 
by-layer deposition on a nanometer scale have attracted much attention due to 
the potential applications of large magnetoresistance. GMR and Tunnel mag- 
netoresistance (TMR) [2-4] discoveries have raised physical issues regarding the 
interplay between the spin and charge of an electron. From the industrial point 
of view, the control of this interplay is important because larger magnetoresis- 
tance is key to higher performance of data storage applications [5-8]. 

Progress in microfabrication technique has enabled us to study characteristic 
phenomena in laterally small sizes, and it has long been known in metallic or 
semiconductor systems, nanometer-sized structures provide insight into meso- 
scopic phenomena such as quantization of conductance, single-electron tunnel- 
ing, discreteness of the chemical potentials, and so on. To study a mesoscopic 
property, we must reduce the size of a structure within a characteristic length ac- 
cording to its feature and temperature. For example, to study Coulomb blockade 
and resulting single-electron tunneling (SET) properties at room temperature, 
the size of a particle should be less than 3 nm (embedded in a vacuum) [9]. 
A microstructure a few nanometer in size, which is under the microfabrication 
limit, is generally prepared by a bottom-up process. 

Returning to the magnetic systems, and in view of GMR or TMR in layered 
nanostructures, it has naturally been anticipated that spin-dependent transport 
in three-dimensional magnetic nanostructures would gives rise to further novel 
phenomena. Our interest emerged especially in the interplay of spin-dependent 
tunneling (SDT) and SET in magnetic nanoparticles. TMR, which was dis- 
covered in 1995[3, 4], has normally been studied in magnetic tunnel junctions 
(MTJs) consisting of a 1-3 nm-thick insulating barrier layer sandwiched between 
upper and lower ferromagnetic electrodes, and its ratio is given by the differ- 
ence in the electrical resistance between parallel and anti-parallel alignments of 
magnetization vectors. The TMR ratio of a macroscopic MTJ with an amor- 
phous insulating barrier is almost consistent with that proposed by Julliere [10] 
in terms of the spin polarizations of the electrodes. In the case that the size of 
the MTJs is reduced to nanometer-scale, the interplay of SDT and SET causes 
a peculiarity in the TMR ratio different from that of the Julliere model. 

Here we report characteristic magnetotranport behaviors arising from the 
interplay of SDT and SET in samples using magnetic nanoparticles. We em- 
ployed an insulating granular film consisting of a few nanometer-sized magnetic 
particles embedded in an insulating matrix. The size of the particles was so 
small that their charging energy exceeds 30 meV (^300 K), in which case SET 
is expected to occur between pairs of neighboring particles at low temperatures. 
In a granular film of macroscopic size containing a large number of particles, 
however, SET phenomena between neighboring particles average out, due to the 
large distribution of particle sizes and interparticle distances. In order to observe 
SET phenomena such as Coulomb staircase and Coulomb blockade with a clear 
threshold, the tunnel paths must be restricted from among the vast number of 
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possible paths by uniting them with microfabricated electrodes. 

First, we describe characteristic magnetotransport behavior in non- microfabricated 
granular films. Although SET properties average out in non-microfabricated 
samples, the higher-order tunneling process in a Coulomb blockade regime, or 
so-called cotunneling, gives rise to enhanced TMR at a low-temperature and 
low-bias-voltage region. Next, we describe magnetotranport behavior for mi- 
crofabricated granular films. The microstructures using granular film we pre- 
pared in this study consisted of either a granular nanobridge or a granular 
CPP (current-perpcndicular-to-plane) structure. In the former, or lateral-type 
microstructure, a granular film was placed in a nanogap of microfabricated elec- 
trodes. We found enhanced TMR associated with Coulomb blockade in this 
structure. The latter was a nanopillar of layered structure: in this case, a gran- 
ular film was sandwiched between electrodes. In the CPP structure, we found 
interesting TMR behavior due to SET and spin accumulation in nanoparticles; 
in addition, it was elucidated that spin relaxation time is much enhanced in the 
Co nanoparticles. 
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2. Basic concepts 



2.1 Spin-dependent tunnel magnetoresistance 

Tunnel magnetoresistance (TMR) originates from spin-dependent tunneling 
in an FM / I / FM magnetic tunnel junction (MTJ), where the FMs are ferro- 
magnetic layers and the I is a tunnel barrier that is a few nm thick. In 1975, 
Julliere formulated a model for the difference in conductance between the par- 
allel and anti-parallel configurations in the two FM electrodes, FM1 and FM2. 
TMR is defined as: 

TMR = ^ = *2^3l = (2.!) 

Up Hp <^AP 

where Rp(ap) is the resistance and Gp(ap) is the conductance for the P (AP) 
alignment [10]. G is given from the sum of the G for each spin channel a, and 
G a is described by 

G ff cx^£ lCTJ D 2(T , (2.2) 

where £ , i(2) (T gs the density of states of the a spin channel at the Fermi level in 
FM1(FM2). The total conductance in the P alignment is given by 

G P = G P + Gp cx D M1 D M 2 + D ml D m2 , (2.3) 

and in the AP alignment 

Gap = G\ P + G l AP cx D M1 D m2 + D rnl D M2 , (2.4) 

where Dui and D m i are the densities of states for the majority and minority 
spin bands in the i-th FM electrode, respectively. The spin polarization of the 
i-th electrode is defined by 

p _ = Dm- D mi _ (2 5) 

Dmi + D m i 

Using this definition, Eq. (2.1) can be expressed as 

TMR = Gp -° AP = 2PlP2 . (2.6) 
Gap I-P1P2 1 7 

According to Eq. (2.6), which was predicted by Julliere, a TMR of about 28 
% is expected assuming both electrodes are Co with P ~ 0.35. This prediction 
of the Julliere model almost exactly reflects the experimental results for MTJs 
with Al-0 amorphous barriers. 

In the case of insulating granular systems, we normally define the TMR 
ratio as TMR = AR/ R max . The TMR ratio should be smaller than those of 
MTJs because of the random orientation of the magnetization (M) at the zero 
field resulting in an incomplete anti-parallel configuration. Inoue and Maekawa 



6 



considered the relative magnetization of a granular system m — M / M s , where 
M s is the saturation magnetization, and derived the TMR as [11]: 

m 2 P 2 

TMR = ' (2.7) 
1 + m z P z 

Assuming Co nanoparticles with P~0.35, the saturated (m = 1) TMR ratio 
P 2 /(l+P 2 ) is evaluated to be about 11 %. 



2.1 Charging effect of nanoparticles 



The charging effect of nanoparticles leads to a Coulomb blockade; then the 
single-electron charge plays a significant role in electron transport [9, 12, 13]. 
A simplified model is shown in Fig. 2-1: an individual particle, isolated in 
an insulating barrier, is located between source and drain electrodes, and the 
tunneling electron travels from the source to the drain through the particle. 
The tunneling of electrons can be inhibited, and the current does not flow at 
small bias voltages if the electrostatic energy e 2 /2C of a single excess electron 
on the island is much larger than the thermal energy ksT, where C is the 
capacitance of the island. The suppression of current at small bias voltages is 
called a Coulomb blockade. When the bias voltage increases and exceeds the 
threshold Vth = e/2C, the current starts to increase. If one junction resistance 
is similar to that of other (i?i w R2), the current increases smoothly with bias 
voltage (see Fig. 2-2 (a)). On the other hand, if the difference between the two 
junction resistances is very large (i?i << i?2or Ri » R 2 ), the current increases 
stepwise with bias voltage depending on the number of electrons accumulated 
on the island (see Fig. 2-2 (b)). The step-like structure in current vs. bias- 
voltage (I — Vb) characteristics is called the Coulomb staircase. The Coulomb 
blockade and the Coulomb staircase are representative phenomena of single- 
electron tunneling (SET). In practice, even for multiple junctions including not 
only one but some islands between electrodes, SET phenomena can be observed. 

For the appearance of SET phenomena, the following two conditions are 
generally essential. First, as mentioned above, the charging energy E c required 
to add an electron to an island must far exceed the thermal energy; i. e., 

E c » k B T. (2.8) 

Second, the junction resistance Rj must exceed the resistance quantum Rk = 
h/e 2 w25.8 kQ; i. e., 

Rj » Rk, (2.9) 

which ensures that the wave function of an excess electron on an island is local- 
ized there. Without this condition, the electron would be delocalized through 
the island, permitting the transport. SET phenomena in nonmagnetic systems, 
where both islands and electrodes are nonmagnetic, have been extensively in- 
vestigated, both theoretically and experimentally [9]. The islands must be small 
enough to satisfy Eq. (2.8), since E c increases as the island size decreases. The 
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E c of an isolated island is described as 

E c = e 2 /4ned, (2.10) 

where d is the diameter assuming the island is a sphere, and e is the dielectric 
constant of an insulator around the island. Although E c is somewhat modified 
by the configuration of surrounding islands and electrodes in a real system, we 
may consider that E c is roughly proportional to 1/d. The E c for submicron- 
sized islands is generally of the order of l(T 4 -l(r 3 cV (1-10 K). Thus, SET 
phenomena can be observed only at very low temperatures. 

There are two distinct tunneling processes by which electrons are transferred 
between the electrodes via a small island: sequential tunneling and cotunneling. 
In sequential tunneling, there is no correlation between tunneling events into 
and out of the island. In this process, electron tunneling in either of the two 
junctions causes an increase of charging energy, and is suppressed by a Coulomb 
blockade. However, cotunneling occurs when tunnel resistance is not so high (1 
or 2 orders higher) compared to the resistance quantum (Rk — 25.8 kO). In 
the cotunneling process, two electrons tunnel in a correlated fashion, i.e., an 
electron tunnels into the island while the second electron simultaneously leaves 
the island through the other junction; the island is only virtually charged in 
this cotunneling process, and therefore there is no increase in charging energy 
in the overall tunneling process. Therefore, the cotunneling process contributes 
to a finite current even in a Coulomb blockade regime. 

2.2 Spin-dependent single-electron tunneling 

In the last section, we neglected the spin degree of freedom. In the mag- 
netic tunnel junctions, the tunnel current depends on the relative orientation of 
magnetizations in the ferromagnetic electrodes. Assuming a double tunnel junc- 
tion model with a ferromagnetic nanoparticle put between ferromagnetic elec- 
trodes, it was expected that the interplay of spin-dependent tunneling (SDT) 
and single-electron tunneling (SET), so called spin-dependent single-electron 
tunneling (SD-SET), would give rise to remarkable tunnel magnetoresistance 
(TMR) phenomena. 

Magnetotransport properties have been studied in both sequential tunneling 
and cotunneling regimes. In the limit of the sequential regime, phenomena are 
discussed in the framework of the orthodox theory of SET. Pioneering theoret- 
ical studies were performed by Barnas and Fert [14, 15], and Majumdar and 
Hershfiled [16] in 1998. They predicted that novel TMR behaviors such as mod- 
ification of the spin-dependent tunneling probability associated with Coulomb 
staircase give rise to a TMR peak around the step point of the staircase. A peak 
appears at each step resulting in the oscillation of TMR as a function of the bias- 
voltage. Following their studies, further analyses were made by several authors 
[17-36]. Some studies have considered SD-SET with additional effects such as 
spin accumulation [17-30, 36] and a discrete energy spectrum [20, 32-34] in the 
island. Characteristic TMR behavior due to the spin accumulation effect was 
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predicted: the difference in the spin-splitting of the chemical potential between 
parallel and anti-parallel alignment of magnetic vectors gives rise to an alter- 
nate sign change of TMR [15, 27, 28, 30, 37]. The spin accumulation effect has 
been mainly considered in an SET device model with a nonmagnetic nano- island 
placed between ferromagnetic electrodes. Although TMR should not appear in 
the framework of the Julliere model, nonequilibrium magnetic polarization due 
to spin accumulation leads to nonzero TMR with periodic oscillation. Studies 
concerning spin accumulation have also examined the magnctotransport prop- 
erties for various spin relaxation times in the island because the accumulation 
occurs when the spin of an electron entering a particle does not flip until the 
next electron arrives. 

On the other hand, in the case of the cotunneling process [38] , magnetotrans- 
port phenomena have been analyzed in a Coulomb blockade region [36, 39-41]. 
In contrast to those in sequential tunneling, which occur outside of a Coulomb 
blockade regime and with in the limit of high tunnel resistance, cotunneling 
gives a dominant contribution when the tunnel resistances are no more than 1 
or 2 orders higher than the resistance quantum. Assuming a double junction, 
the total tunnel resistance for sequential tunneling is proportional to the sum 
of the resistances, while that for cotunneling is proportional to their product. 
It has been predicted that this discrepancy would lead to enhanced TMR in a 
cotunneling regime [41]. 

Experimentally, there have been several studies [42-72]. The first experi- 
mental evidence for characteristic SD-SET phenomena was observed by Ono et 
al. in 1995 [42]. They fabricated Ni / NiO / Co / NiO / Co double junctions 
with small contact area (^200nm 2 ) and measured magnetotransport properties 
by changing the gate voltage at a low temperature, below 1 K. In the Coulomb 
blockade region, they found enhanced TMR execededing 40 % in the off-state 
but only reacing 4 % in the on-state. The enhancement was theoretically ex- 
plained by a strong tunneling model in which the tunnel resistances were near 
the resistance quantum [39, 73]. Although the authors of these studies suc- 
cessfully observed enhanced TMR, the operating temperature was lower than 
100 mK because the size of the microfabricated island was larger than 100 nm. 
Schelp et al also found a Coulomb blockade effect on TMR in a sample with 
nanometric Co clusters [55]. They prepared a layered Co / Al-0 / Co clusters 
/ Al-0 / Co sample and observed twice-larger TMR at 4.2 K than that at RT. 
The origin for enhanced TMR was suggested to be the effect from the Coulomb 
blockade in Co clusters. Their observation of the Coulomb blockade and TMR 
even in the sample with a very large contact area was of significance, because 
they showed the potential of studying SD-SET in a layered structure. Follow- 
ing their study, several other SD-SET studies in layered structures have been 
performed [48, 53, 56, 58, 63, 65]. In order to understand the mechanism of 
SD-SET, as mentioned above, it is necessary to divide the properties of SD- 
SET into two categories according to transport mechanisms: cotunneling and 
sequential tunneling. Although there have been some reports on enhanced TMR 
ascribed to finite spin-dependent transport due to the cotunneling process, no 
experimental evidence in a sequential tunneling regime has been reported. We 
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have fabricated appropriate sample structures using Co-Al-0 insulating granu- 
lar films, and elucidated SD-SET phenomena in sequential tunneling [60, 63-65] 
and co-tunneling regimes [62, 66] individually. Our study is the first to reveal 
an oscillatory TMR behavior in a sequential tunneling regime. The details are 
described later. 



2.3 Spin accumulation in nanoparticles 



Spin injection and accumulation were first studied in a layered sample where 
a normal metal (NM) layer was sandwiched between two ferromagnetic (FM) 
layers. A non-equilibrium effect of spin accumulation is generated by a spin- 
polarized current from one FM layer into the NM layer. When the NM layer 
is thinner than the spin-diffusion length, a spin injection signal attributed to 
spin accumulation in the NM layer is detected in the other FM layer according 
to the relative orientation of the two FM layers [74]. Spin accumulation has 
also been studied in sub- micrometer-sized lateral structures [75, 76]. In this 
manuscript, we discuss the spin accumulation effect in an isolated nanoparticlc 
in which spin-polarized current is injected through a tunnel barrier. In such a 
system, spin accumulation in the nanoparticle occurs when the spin relaxation 
time overcomes the interval of successive electron tunneling [15, 19, 21, 23, 25, 
36, 37, 77, 78]. According to the sample structure shown in the Ch. 6, we 
introduce the FM / I / FM nanoparticle / I / NM model (I: insulating barrier) 
as illustrated in Fig. 2-3. As mentioned in Ch. 2, Scc.l, a tunnel conductance 
G a for a spin states ergs expressed by the multiplication of the initial and final 
density of states (DOS) of the a spin band at Fermi level (Ef)- There should 
be a discrepancy between incoming and outgoing numbers of spins at a fixed 
spin state because the DOS at Ep of FM is spin-polarized while that of NM is 
not. This gives rise to the increase (or decrease) of spin population leading to 
the shift of the chemical potential AE F for a a spin state. In order to maintain 
charge neutrality in the particle, the chemical potential of another a shifts in 
the opposite direction, satisfying D^AE F = —D^AEp. Although the degree 
of the shift depends on the net DOS of the particle, AE F is not the same as 
AEp because the net DOS is spin-polarized around E F . In order to maintain 
spin-conservation conditions, the following expression is applied: 

(/lW2 " } = ^AE F (2.11) 

e tsf 

where Jj l(7 is the current at the i th junction (i= 1, 2) for spin a, O is the volume 
of the particle, tsf is the mean spin relaxation time in a nanoparticle defined as 
Tg F = (t^ 1 +t j ^ 1 )/2, and To- is the spin relaxation time of electrons with spine. 
When we consider the spin accumulation effect in SET regime, the bias-voltage 
dependence of Ai?£.shows a sawlike oscillation with a period close to that of the 
Coulomb staircase. In the case of the Co nanoparticle, characteristic oscillation 
gives rise to modifications in the current - bias voltage curve resulting in novel 
TMR behavior. Details will be shown in Ch. 6. 
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3. Insulating granular film 

3.1 Materials, sample preparation and microstructure 

Insulating granular films consist of small metallic particles and an insulating 
matrix. If the composition of the metal is lower than the percolation limit, 
the transport is dominated by the tunneling of electrons between particles. If 
the particles are magnetic, TMR arises because the magnetization vectors on 
particles, which are not aligned at low applied fields, become more aligned as 
the applied field increases, leading to a decrease in resistivity. Pioneering works 
on TMR in insulating granular films were published by Gittleman et al. .[79] 
and by Helman and Abeles [80], both of whom used Ni-Si-0 films. However, the 
magnitudes of TMR were very small. In 1994 Fujimori et al. used a Co-Al-0 
granular film and reported a large TMR, which reached 10 % even at rooom 
temperature [81]. This ratio was found to increase with decreasing temperature, 
exceeding 20 % at low temperatures [62, 82]. 

Insulating granular films are easily prepared by sputtering and evaporation 
techniques. As shown in Fig. 3-1, the primary methods we used included (a) 
reactive sputtering, (b) sputtering with a composition target, and (c) tandem 
deposition with plural targets. We show the details of these methods for the 
preparation of a Co-Al-0 granular film. In the case of (a), a Co-Al alloy target 
was sputtered in an Ar+02 atmosphere; then Al was selectively oxidized. In 
the case of (b), AI2O3 sheets are were placed over a Co target, and sputtering 
was performed in an Ar atmosphere. In the case of (c), there were two cathodes 
with targets, Co and AI2O3, and the rotation of the substrate holder caused 
the alternation of their deposition. In all cases, a granular structure comprising 
nanometer-sized Co metal particles embedded in an insulating Al-0 matrix is 
formed on substrates. 

In this manuscript, we employed Co-Al-0 films prepared by the reactive 
sputtering method. Hereafter the Co-Al-0 granular film containing x at.% Co, 
y at% Al and z at.% O is denoted Co^Al^O^. The compositions are deter- 
mined by Rutherford backscattering spectroscopy (RBS) analysis. The atomic 
composition of Co was roughly controlled by that in the Co-Al target. On the 
other hand, the composition of the film, at its fixed target composition, is rel- 
atively insensitive to the oxygen gas flow ratio. For example, in the case of the 
C036AI22O42 film, the target composition was fixed to be C025AI75 while the 
oxygen gas flow ratio was allowed to range from 1.6 % to 3.0 % (at total gas 
pressure ~ lmTorr). 

Figs. 3-2 show (a) plan view and (b) cross-sectional transmission electron 
microscopy (TEM) micrographs for Co46Al 19 035 film [83]. The film has an 
isotropic granular structure consisting of Co particles 2-3 nm in diameter (dark 
spheres) and intergranular Al-0 about 1 nm in thickness (white channels). Fig. 
3-3 shows a high resolution TEM micrograph for C052AI20O28 film, indicating 
that a crystalline Co particle is surrounded with non-stoichiometric aluminum 
oxide with amorphous structure. 
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3.2 Electric transport and magnetic properties in insulating gran- 
ular systems 

As mentioned above, the electric charging effect of small islands causes a 
Coulomb blockade when the charging energy overcomes the thermal and/or 
bias potentials. An insulating granular system consists of metallic nanoparti- 
cles embedded in an insulating matrix. Each nanoparticle is electrically isolated; 
their average size is a few nanometers while their charging energy exceeds 100 
K. The electric transport in insulating granular systems represents a character- 
istic dependence on temperature. Fig. 3-4 shows the temperature dependence of 
electrical resistivitygp for Co-Al-0 films of different compositions: C054AI21O25, 
C052AI20O28, C046AI19O35, and C036AI22O42. The plot of In p versus T -1 / 2 is 
approximately linear for all the films. The relationship of In/? cx l/v^was first 
derived by Sheng et al. [84], considering electron tunneling, charging effect and 
particle size distribution. They analyzed its temperature dependence assuming 
a model in which particles on each contribution path have the same size d and 
are separated by a barrier thickness s, keeping their ratio s/d (or equivalently 
E c s) constant for a given composition. In this model, the electrical conductiv- 
ity is dominated by tunneling between small particles at high temperature and 
between large ones at low temperatures because larger ones cause a Coulomb 
blockade as the temperature decrases. 

For a conventional sandwich-type tunnel junction, the tunnel conductance is 
temperature independent and proportional to exp{- 2(27r/ft,)(2?7j0) 1 / 2 s}, where 
h is the Planck's constant, m the effective electron mass, (f> the effective barrier 
height, and s the barrier width. In the case of an insulating granular system, the 
existence of finite charging energy (E c ) gives rise to the particular temperature 
dependence of p. The number density of the charge carrier, whose generation 
requires an E c , is proportional to the Boltzmann factor expj.Ec/2/cB7 1 }, where 
E c is the Coulomb energy, ks the Boltzmann's constant, and T the absolute 
temperature. The final formula about p was obtained as 

p = Po {cxp{2(C/A :B ) 1 / 2 r- 1 / 2 }, (3.1) 

where C — (2tt /h){2m4>) 1 ^ 2 sE c which was named the tunnel activation energy. 
Eq. (3-1) can be transformed to 

lnp = 2(C/fc i ?) 1/2 T- 1 / 2 + const. (3.2) 

The linear relation found experimentally in Fig. 3-4 agrees with Eq. (3-2). 
From the gradient of the log p versus T -1 / 2 , C is estimated to be 30 meV for 
the Co 36 Al22 42film. 

Figs. 3-5 (a) and (b) show M — H and M — T curves, respectively, for the 
Co 36 Al 2 2042 film. The M-H curves at RT show no remanence at zero field and 
unsaturated behavior in high fields. Large thermal hysteresis is observed be- 
tween zero-field cool (ZFC) and field cool (FC) curves, and it is suggested that 
particles are superparamagnetic at room temperature and that thermal fluctu- 
ation of the magnetic moments is blocked at 4.5 K. Figs. 3-6 (a) and (b) show 



12 



MR curves for C036AI22O42 measured at RT and 4.2 K, respectively. Large 
values of MR exceeding 10 % are observed. At RT MR does not completely 
saturate even at 80 kOe; on the other hand, MR nearly saturates above 20 kOc 
at 4.2 K. The coercivity and hysteresis of the MR curve appearing at 4.2 K 
corresponds to the magnetization curve shown in Fig. 3-5 (a). 

3.3 Particle size distribution in insulating granular systems 



The sizes of nanoparticles in insulating granular films are distributed; in- 
formation on the particle size distribution is essential to understanding the 
mechanism of TMR phenomena [85]. Ohnuma ct al. [83, 86] reported the mi- 
crostrustures in Co-Al-0 granular films investigated by high-resolution trans- 
mission electron microscopy (HRTEM) to obtain the distributions of particle 
sizes and interparticle distances. We have evaluated particle size distributions in 
C036AI22O42 granular films by fitting the magnetization curves to the Langevin 
function which describes superparamagnetic behavior. 

Magnetization in the superparamagnetic state can be described by the Langevin 
function. If we assume that all the Co particles have a spherical body with the 
same diameter d and if the anisotropy energy is negligible, the magnetization 
M is described by 

M 1 

= L (a) = coth a (3.3) 



M s 



a 



with 



a = 



HH _ M Co H 4tt 
k B T ~ k B T 3~ 



(3.4) 



where Ms is the saturation magnetization of the sample, MCo is the magne- 
tization of Co particles and H is the applied field. In real systems, d has a 
distribution; we use a log-normal distribution function (LNDF) in d. LNDF is 
described by 



f(d) 



1 



2tt ln< 



■ exp 



(hid — lad m y 
21n 2 CT 



(3.5) 



where the parameters dm and a represent the statistical median and the ge- 
ometric standard deviation, respectively. The Langevin function taking the 
log-normal particle size distribution into account is described by 



Mcof (f ) /K-)Mn,-) 



(3.6) 



with 



M Co H An 
k B T ~~3 



(3.7) 



Fig. 3-7 (a) shows the calculated magnetization curves compared to the experi- 
ment for a C036AI22O42 film at T = 200 and 293 K; the particle size distribution 
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which gives the best fit is shown in Fig. 3-7 (b). The sizes of particles are dis- 
tributed mostly in the range of 1.0 to 5.0 nm in diameter and dmag is estimated 
to be 2.45 nm. This dmag is consistent with that evaluated from the TEM 
image (Item (= 2.20 nm). 
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4. Spin-dependent higher order tunneling in Coulomb blockade 
regime 

In this chapter, we report on the temperature and bias- voltage dependence of 
MR in insulating Co-Al-0 granular films without performing microfabrication. 
The MR exhibits an anomalous increase at low temperatures but no significant 
indication of change with bias voltage. We show that the anomalous increase 
of the MR indicates evidence for higher-order tunneling (cotunneling) between 
large granules through intervening small granules [36, 62]. We emphasize that 
the existence of higher-order tunneling is a natural consequence of the granular 
structure, since a broad distribution of granule size is an intrinsic property of 
granular systems. 

4.1 Results 

Fig. 4-1 shows the temperature dependence of MR for Co-Al-0 films. It 
is clearly seen that the MR is remarkably enhanced at low temperatures while 
it is nearly independent of temperature above ~100 K. For C036AI22O42, the 
MR below 3 K is anomalously large and reaches more than twice the value 
given by PcolO-^^Co) [H]> where the formula is half of that for magnetic tun- 
nel junctions (MTJ) because of the difference between random and antiparallel 
alignment of magnetic moments. In the case of MTJ with an Al-0 barrier, 
the temperature dependence of MR is discussed on the basis of magnetic im- 
purity or magnon scattering. However, it is considered that magnetic impurity 
or magnon scattering does not give rise to the plateau in the temperature de- 
pendence of MR, as in the present result observed above ~100 K. Helman and 
Abeles [80] proposed a theory of spin-dependent tunneling in insulating granular 
systems and predicted the temperature dependence 1/T for MR. However, the 
dependence 1/T does not fit the present results. 

Transport properties were measured in the current-perpendicular-to-plane 
(CPP) geometry as shown in Figs. 4-2 where a 1-^.m-thick C036AI22O42 gran- 
ular film was sandwiched between upper and lower Au-Cr electrodes. Figs. 4-3 
(a) and (b) show p and MR at 4.2 K, respectively, as functions of bias voltage 
Vb for a C036AI22O42 film, p decreased rapidly by 3 orders of magnitude as 
the bias voltage increased from Vb = up to 600 mV. Nevertheless, the mag- 
nitude of the enhanced MR was almost constant. This is in clear contrast with 
the case of MTJ of macroscopic size, where both MR and p decrease gradually 
with increasing bias voltage. Furthermore, the bias-voltage dependence of MR 
is much smaller than the temperature dependence of MR in Fig. 4-1. We can 
consider that about 200-300 Co granules exist in the direction normal to the 
film plane between the upper and the lower electrodes, assuming an average 
particle size (d) of 2.5 nm and interparticle distance (s) of 1 nm. Therefore, 
the applied bias voltage per one microjunction consisting of two neighboring Co 
granules may be estimated to be 2-3 mV at Vb = 600 mV, which corresponds 
to 20-30 K in temperature. As seen in Fig. 4-1, the enhanced MR decreases 
rapidly with increasing temperature and becomes flat around 20-30 K, while it 
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is independent of Vb at least up to 600 mV (2-3 mV at neighboring ones). 



4.2 Theoretical explanation 

In granular systems with a broad distribution in particle size [85] , it is highly 
probable that large particles are well separated from each other due to their 
low number density (i.e., the larger the granule size, the more separated the 
granules), and there may be a number of smaller granules between large granules 
as shown in Fig. 4-4 (a). To model the structural feature of granular systems 
we assume that large granules with size n(d) and charging energy (E c )/n are 
separated by an array of n particles with average size 

(d) 

and charging energy (E c ) on a conduction path, as shown in Fig. 4-4 (b). 
4.2.1 Temperature dependence of MR 

We first calculate the temperature dependence of the conductivity <r(T) at 
zero- bias voltage (Vb = 0). The tunnel current at the zero-bias limit is dominated 
by thermally activated charge carriers. In the case of the conduction path in 
Fig. 4(b), the carriers mostly occupy the large particle of charging energy 
(E c }/n in a probability proportional to the Boltzmann factor exp[— (E c ) /2nT] 
in units of ks = 1- Since the large particles are separated by the smaller 
ones, the ordinary tunneling of an electron from a large particle to a small one 
increases the charging energy by 8E C <~ I/2 (E c ) /(l + 1/n); thus is suppressed 
by the Coulomb blockade at low temperatures T < SE C . In this regime, 
the dominant contribution toger(T) comes from higher-order processes of spin- 
dependent tunneling where the carrier is transferred from the charged large 
particle to the neighboring neutral large particle through an array of small 
particles, using successive tunneling of single electrons, i.e., the cotunncling of 
(n+1) electrons. Figure 4-4 (a) shows an example of the third-order process 
(n=2). Summing up all of these higher-order processes, we have 

a(T) OC Y e -<£c>/2nT [(1 + p2 m 2 )e -2 K s' ]n+ l I W^f 

„ \(SE C ) 2 +^(T) 

Here, [ — ] is the spin-dependent tunneling probability between the neighboring 
particles, m = M/M s is the magnetization normalized to the saturation mag- 
netization M s , k is the tunneling parameter related to the barrier height, and 
s' = 2n (s) /(n + 1) with (s) being the mean separation of particles with size(d). 
The factor ( — )™ represents the finite temperature effect by which electrons (or 
holes) in the energy interval of ttT around the Fermi level participate in the 
intermediate states of the higher-order process [87], and j(T) is the decay rate 
given by "f(T) w gT with g being a constant .[88]. The function f(n) represents 
the distribution of the conduction paths. In Eq. (4.1), exp [Ann (s) — (E c ) /2nT] 

1/2 

is a peaked function of n and has its maximum at n* = ((E c ) / 8R (s) T) 1 , where 



f(n) (4.1) 
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k/k w 1 + (1/4/t (s)) In [(.g/7r) 2 + ((£ c ) /2ttT) 2 ] . The existence of higher-order 
tunneling processes with different orders is a natural consequence of the granu- 
lar structure since a broad distribution of granule sizes is an intrinsic property 
of granular systems. Namely, cotunncling processes with n = 1 occur in some 
places, while those with n — 2 occur in other places, thus n (and consequently 
n*) can be treated as a continuous variable at low temperatures (T « (E c )). 
Replacing the summation by the integration in Eq. (4.1) and using the method 
of steepest descent [89], we obtain 



a (T) oc (1 + P 2 m 2 ) 



2\ n* + l 



k{s) 



f(n*) exp 



2k (s) (E c ) 



(4.2) 



(4.3) 



In Fig. 4-5, the calculated resistivity for in = is shown by the solid lines. Here 
and hereafter, we assume / (n*) ex 1/n*, and take 2k (s) = 3, g = 0.3, and the 
values of (E c ) are estimated to be 9 K for C054AI21O25, 18 K for C052AI20O28, 
25 K for C046AI19O35, and 110 K for C036AI22O42. 

Because of the higher-order processes, the spin-dependent part of cr(T) in 
Eq. (4.2) is amplified to the (n* + l)th power of (1 + P 2 m 2 ), so that a(T) is 
sensitive to the applied magnetic field since m varies from m = to m = 1 (the 
fully magnetized state) by application of the magnetic field. Using Eq. (2) the 
MR, Ap/po = 1 - [a (T)] m=0 /a (T), is expressed as 

^ = !-(! + m 2 P 2 )- ( "* +1) . 
Po 

The calculated MR is shown by the solid curves in Fig. 4-1, where the value of 
P is chosen to fit the experimental data. For small P 2 , Eq. (3) is approximated 
to be 

with C = (E c ) /8k (s) being constant. Eq. (4.4) indicates an anomalous in- 
crease of Ap/po at low temperatures due to the onset of higher-order processes 
between larger granules, i.e., n* oc 1/VT. At T = 2 K, n* takes the value of 1.6, 
so that one or two small granules intervene between larger ones in the higher- 
order processes. As seen in curve a, the MR grows rapidly around 10 K well 
below E c = 110 K. Similar behavior is seen in a double-junction system [41] [36]. 



(4.4) 



4.2.2 Bias voltage dependence of MR 

We next calculate the bias-voltage dependence of conductivity a (14) m the 
Coulomb blockade regime. When a finite voltage 14 is applied to the granular 
system, the voltage drop A14 between the large granules in the model system 
of Fig. 4-4 (b) is given by AV4 = (2n/N g ) 14, where N g is the average number 
of glanules along a conduction path. cr(Vb) at finite temperatures, the factor 

(ttT) 2 ™ in Eq. (4.1), is replaced by (ttT) 2 + {2eV b /N g ) 2 [87]. We obtain 



er(14) oc J^V 



(E c )/2nT 



[(I + P 2 m 2 )e- 2KS } 



2ns f ]n-\-l 



2eV b 
N g nT 



f(n). (4.5) 
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Following the same procedure as in deriving a(T) in Eq. (4.2), we obtain the 
bias -dependence of the conductivity 



The a(Vb) exhibits a power low dependence (1/VJ,) 2 "* for T < eV h jN g . 

In Fig. 4-2 (a), we show the calculated resistivity po (Vb) — l/cr (Vb) by the 
solid curves for T = 4.2 K and N g — 140. The steep decrease of the calculated 
resistivity is in good agreement with that of the experimental data. In Fig. 
4-3 (b), the calculated MR is shown by the solid curve. The enhanced MR is 
maintained upon application of higher voltages, which is consistent with the 
experimental result. The constant MR may originate from the large number of 
granules along the conduction paths in the granular films, in which the voltage 
drop between neighboring granules ~ V/N g is small for a large value of N g . V / 
N g at Vb = 500 mV is about 3 mV for N g = 140 and its corresponding temper- 
ature is ~ 40 K, which is smaller than that of the charging energy (E c ) ~ 110 
K. We note that the enhanced MR is nearly constant up to 500 meV, whereas 
the corresponding resistance is reduced by several orders of magnitude. This is 
in contrast with ferromagnetic tunnel junctions of macroscopic size, where both 
the MR and the resistance decrease gradually with increasing bias voltage. 




(4.6) 
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5. Spin-dependent single-electron tunneling in microfabricated 
structures 

5.1 Current path restriction to observe spin-dependent single-electron 
tunneling in sequential tunneling regime 

Insulating granular films consisting of nanometer-sized magnetic metallic 
particles embedded in an insulating matrix are useful for the study of spin- 
dependent SET (SD-SET) phenomena. The properties of SD-SET are roughly 
divided into two categories according to transport mechanisms: cotunneling 
and sequential tunneling. In the last chapter, we discussed one characteristic 
SD-SET behavior due to cotunneling in the Coulomb blockade regime. Here 
we discuss about the SD-SET behaviors in a sequential tunneling regime: in 
other words, TMR behaviors associated with a Coulomb staircase and/or clear 
Coulomb threshold, in microfabricated samples. 

In granular film of a macroscopic size containing a large number of particles, 
however, SET phenomena, represented by Coulomb staircase and so on, are 
averaged out due to the large distributions of particle sizes and interparticle 
distances [62, 83, 85]. The tunnel paths should be restricted in order to clearly 
observe the SET phenomena. A simple method to restrict the tunnel paths is 
to use scanning tunneling microscopy (STM). The tunnel path on the surface is 
limited to only one particle just below the STM tip. We observed clear Coulomb 
staircases in the I — Vb measurements for Co-Al-0 granular films even at room 
temperature [90, 91]. Figs. 5-1 (a) and (b) show typical examples of an STM 
topographic image and an I — Vb curve, respectively, for a C036AI22O42 film. 

A more advantageous method than STM for a variety of measurements 
and applications is to fabricate a device structure consisting of a small part 
of granular film with nanofabricated electrodes. In this study we have fab- 
ricated two types of device structures with Co-Al-0 granular films using fo- 
cused ion-beam (FIB) milling or electron-beam lithography techniques. One 
is a granular nanobridge structure [60, 67]: point-shaped electrodes separated 
by a very narrow lateral gap filled with Co-Al-0 granular film. The other is a 
current-perpendicular-to-plane (CPP) geometry structure [63-65]: a thin Co-Al- 
O granular film sandwiched by ferromagnetic electrodes with the current flowing 
in the direction perpendicular to the film plane through a few Co particles. We 
measured the I — Vb curves in these samples, and found the enhancement and 
oscillation of TMR due to spin-dependent SET in a sequential tunneling regime. 

5.2 Results in nanobridge structures 
5.2.1 Experimental procedures and results 

In this section, we report enhanced TMR just above the Coulomb blockade 
threshold in a sequential tunneling regime [60]. We fabricated point-shaped 
electrodes separated by a very narrow lateral gap filled with insulating granular 
film; we call the resulting structures "granular nanobridges" . As mentioned in 
the last section, the mechanism for TMR enhancement is different from that 
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for higher-order tunneling, because this enhancement in nanobridgc samples is 
caused outside the Coulomb blockade regime. We apply the orthodox theory of 
SET [9] and explain that enhanced TMR is brought about by the modification 
of the detailed balance of particle charges by the external magnetic field [14, 16, 
91]. 

A schematic view of a typical sample is shown in Fig. 5-2 (a). An insulating 
granular nanobridge was fabricated on a Corning no. 7059 glass substrate as 
follows: a 15 nm-thick-NbZrSi amorphous layer was deposited by rf sputtering, 
and was formed into source and drain electrodes by FIB etching using 30 kV 
gallium ions (Seiko Instruments Inc., SMI 9200). The electrodes separated by a 
gap with a length (I), i.e., gap separation, of 30 nm and a width (w) of 60 nm are 
shown in Fig. 5-2 (b). Deep trenches (60 nm wide and 200 nm deep) were formed 
beside the gap by FIB etching to avoid the formation of unnecessary current 
paths outside the gap. A 7.5-nm-thick Co-Al-0 granular film was deposited on 
the patterned surface by reactive rf sputtering, and the gap was filled by the Co- 
Al-0 film. The aspect ratio of trenches was so high that the trenches were not 
filled with the Co-Al-0 film. The composition of Co-Al-0 was determined to be 
C036AI22O42 by RBS analysis. The average size of Co particles was estimated 
to be about 2.5 nm from the analysis of the superparamagnetic behavior and 
TEM observation [85]. The characteristic sizes of granular nanobridges, i.e., 
w, I, and thickness (i), varied in the range of 60 - 700 nm, 30 - 70 nm and 
5-30 nm, respectively. I — Vb characteristics were measured at 4.2 K using 
an electrometer (Keithley 6514) with a two-terminal arrangement. TMR (= 
AR/R max ) was evaluated from the difference between the I —Vb curves at the 
applied field H = and 10 kOe. 

Fig. 5-3 (a) shows the I — Vb curves at H=0 (solid lines) and H=10 kOe 
(dashed lines) for the sample with w—60 nm, Z=30 nm and t=7.5 nm. Here, 
the threshold voltage Vth was defined as that below which the current was zero 
within an accuracy of 100 fA, and it was approximately 1.5 V in this case. In 
the range of |T4| < Vh, a Coulomb blockade occurred, and the current increased 
rapidly when \V,\ exceeded V t h- It is noted that for Co-Al-0 granular films of 
macroscopic sizes, the Coulomb blockade has not clearly been observed, because 
the macroscopic sample contains a lot of Co particles with a broad distribution 
of sizes and the tunneling of electrons between large particles with small E c 
could start at small voltages. In the granular nanobridge, however, the tunnel 
paths are so limited that the Coulomb blockade is remarkable. 

Fig. 5-3 (b) shows the Vb dependence of TMR. TMR depends strongly 
on Vb- For \Vb\ <4.0 V, TMR increases with decreasing |Vj,| and reaches a 
maximum value larger than 30 % at the voltage slightly above V t h (~1.5 V). 
In the Coulomb blockade region, i.e., \V,\ < V t h (hatching area), there is little 
quantitative reliability of the measurements because the current is very low 
(<100 fA). For \V b \ > 4.0 V, on the other hand, TMR is about 8 % showing no 
large change with Vb, although weak oscillatory behavior seems to exist in the 
V b dependence of TMR. 

Similar results have been obtained in other granular nanobridges of different 
sizes. Figs. 5-4 (a) and (b) show the I — Vb curve and the Vb dependence of 
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TMR, respectively, for the sample with w=700 nm, 1=40 nm and £=15 nm. Vth 
is observed to be 0.5 V, which is lower than that in the sample shown in Fig. 
5-3. Vth shows a tendency to increase as the size of the granular nanobridge 
decreases. The voltage where the TMR shows a maximum, V p , is slightly larger 
than Vth- Fig. 5-5 shows V p vs. V t h in granular nanobridges of different sizes. 
A clear correlation between V p and V t h is seen, suggesting that the enhanced 
TMR is caused by the Coulomb blockade. 

5.2.2 Theoretical explanation 

The orthodox theory of SET [9] can be used to explain the experimental 
results, particularly the enhanced TMR near Vth- In the orthodox theory, the 
tunnel path is modeled as an equivalent classical electrical circuit. For the 
granular nanobridge, we assume a parallel circuit of triple tunnel junctions as 
shown in Fig. 5-6 (a). This is the simplest model to explain our experimental 
results because we need at least two magnetic particles in each series of junc- 
tions to study the spin-dependent transport in nanobridges with nonmagnetic 
electrodes. We neglect the higher-order tunneling process, i.e., cotunneling [41, 
62], because the tunnel resistances between particles and between an electrode 
and a particle are estimated to be about 10 5 times larger than Rq «25.8 kO. 
In order to obtain a stable tunneling current, we constructed a detailed balance 
equation for the probability of states p({rii} a ), which is given in the matrix 
form by 

p = Mp = 0, (5.1) 

where p = (. . . ,p({rii} a ), . . .) T , and M is the transition matrix in the configu- 
ration space constructed by {rii} a with the index a labeling the different charge 
states. The tunneling current through the k th junction is given by 

h=eJ2 a P({n i }a)[rt({n i } a )-r^({n i } a )], (5.2) 

where r^~^ ({rn} a ) cx 1/Rk is the forward (backward) tunneling rate through 
the k th junction with the initial charge state {rii} a . Current conservation re- 
quires that the tunneling current is the same for all the junctions in each 
series. Let us evaluate the tunneling current at the junction between a nonmag- 
netic electrode and a particle, where the tunneling rate \{m} a ) is inde- 
pendent of the magnetic field. The magnetic field dependence of the tunneling 
current comes from the probability p({rii} a ) of the charge state {rn} a , which is 
determined by Eq.(5.2). Since the transition matrix M contains the tunneling 
rates between magnetic particles, the matrix M and therefore the probability 
p({ni} a ) can be modified by applying the magnetic field. For V, just above 
V t h, w e have a few charge states contributing to the tunneling current, and the 
tunneling rates at these charge states differ greatly from each other according 
to thier charging energies. Therefore, the strong modification of the probabil- 
ity p({rii} a ) is made to satisfy the detailed balance equation, and the TMR is 
strongly enhanced just above Vth- This kind of TMR enhancement in double 
tunnel junctions has been studied by Barnas and Fert [14] and by Majumdar 
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and Hcrshficld [16] . They also predicted the oscillating behavior of TMR against 
Vb associated with the Coulomb staircase. No Coulomb staircase appears in the 
I — Vj, curves of granular nanobridges; however, a small-magnitude oscillation of 
TMR could be observed (see Figs. 5-3 and 5-4). There is no strong asymmetry 
in tunnel resistances, i.e., no bottleneck in tunnel paths of granular nanobridges, 
leading to no Coulomb staircase and only weak oscillation of TMR. Moreover, 
we have many junction arrays in a granular nanobridge, and the randomness of 
junction capacitances also muddies the TMR oscillation. 

In order to explain the experimental results for the sample with w=Q0 nm, 
1=30 nm and £=7.5 nm (Fig. 5-3), we considered a parallel circuit of 20 triple 
tunnel junctions and assumed that the tunnel resistance between an electrode 
and a particle is expressed as R ep = (1 ± S)R ep , where S is the deviation from 
the typical value R ep - Other junction parameters such as tunnel resistances 
between particles R PP , junction capacitances C ep , and C pp were also assumed 
to be distributed around the mean values, i.e., R pp = (1 ± S)R PP , C ep = (1 ± 
5)C ep , and C pp = (1 ± S)C' pp . The deviation 5 for each junction parameter was 
randomly chosen within the range of -0.1<<5<0.1. The temperature was set to 
be 4.2 K and the typical value of tunnel resistances for the parallel alignment 
of magnetizations was taken to be 

Rpp = Rep/2- (5-3) 

The tunnel resistance between particles for the antiparallel alignment of magne- 
tizations was larger than that for the parallel alignment and is expressed using 
the spin polarization P as 

R pp = (R ep /2)-(l + P 2 )/(l-P 2 ), (5.4) 

where P is assumed to be 0.42 for Co [92]. The typical values of junction capac- 
itances are taken to be C ep = O.laF and C pp = 0.05aF, which are reasonable 
values considering the average particle size and interparticle distance in Co-Al-0 
granular films [83, 85]. 

The Vb dependence of TMR obtained by the numerical calculation is shown 
in Fig. 5-6 (b). One can see that the theoretical result is in good agreement 
with the experimental one. The TMR is enhanced just above Vth (~1.5 V) 
and decreases with Vb- The randomness of junction capacitances obscures the 
oscillation of the total TMR as shown in Fig. 5-6 (b). We may consider that 
the quantitative difference between Vth an d V p is not essential, but it is caused 
by effects such as the leakage of current through the glass substrate. 

5.3 Results in CPP geometry structures 

In the last section, we fabricated nanobridge structures for combining in- 
sulating granular films with microfabricated electrodes and successfully found 
enhanced TMR at the Coulomb threshold voltage. Proper limitation of the 
number of current paths made it possible to observe spin-dependent SET. How- 
ever, no Coulomb staircase was observed. In an assembly of nanoparticles such 
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as granular films, the Coulomb staircase is expected to appear when the tunnel 
resistance between two neighboring particles or between a particle and an elec- 
trode is much larger than the other resistances in the current path. In other 
words, a bottleneck of tunnel conductance must exist somewhere in the current 
path [9, 36]. In this section, in order to investigate the relationship between 
the Coulomb staircase and TMR, we employed CPP (Current-Perpendicular- 
to-Plane) geometry measurements in Co-Al-0 granular films. A bottlcncch is 
easily added to samples for CPP geometry measurements for the sake of ob- 
serving Coulomb staircases. We prepared CPP geometry samples, where a thin 
Co-Al-0 granular film was sandwiched by ferromagnetic electrodes and a very 
thin Al-0 layer was inserted between the bottom electrode and Co-Al-0 as a 
bottleneck, and measured the current (J) - bias voltage (V&) characteristics. We 
succeeded in observing clear Coulomb staircases due to the current confinement 
among the vast number of channels between the upper and lower electrodes [65] . 

We fabricated CPP geometry samples by a focused ion-beam (FIB) etching 
process. Fig. 5-7 (a) represents a schematic illustration of a sample. Samples 
were prepared on Si/SiC>2 substrates by rf sputtering. A bottom electrode was 
first deposited, and then a thick Al-0 film (40 nm) was deposited onto the 
bottom electrode as an insulating layer using an AI2O3 target. Next, a small 
contact area was made by FIB milling (Fig. 5-7 (b)). A contact area of about 
0.5 x 0.5 /1m 2 was estimated from the scanning ion microscopy image. The 
milling process was carefully performed to leave a very thin Al-0 layer, which 
contributes to forming a bottleneck. After making the contact window, further- 
more, a 1~2 nm thick Al-0 layer was deposited. Consequently, the bottleneck 
is given by the combination of the residual Al-0 and the deposited Al-0 layer. 
A 7^12 nm-thick Co-Al-0 granular layer followed by a top electrode was fi- 
nally deposited. The deposition of Co-Al-0 granular films was done through 
the use of reactive sputtering in a mixture of Ar+02 atmosphere. I — Vb char- 
acteristics were measured at 4.2 K using an electrometer (Keithley 6514) with 
a two-terminal arrangement. TMR (= AR/R max ) was evaluated from the dif- 
ference between the I —Vj, curves at the applied field H = and 10 kOe. The 
external magnetic field was usually applied in the direction parallel to the plane. 

Fig. 5-8 (a) shows I — Vj, curves at H — (solid lines) and H = 10 kOe 
(dashed lines) at 4.2 K. Clear Coulomb staircases are observed for both I —Vt, 
curves. The first three steps from zero bias appear every 20 mV; however, the 
steps at higher bias voltages do not maintain a regular period. In spite of more 
than 10 4 parallel current paths between the electrodes in the contact area (0.5 
x 0.5 /im 2 ), which is much larger than the average Co particle size (2.5 nm), 
clear Coulomb staircases were observed. This suggests that the current at low 
bias voltages preferentially flow along certain restricted local paths where the 
total charging energy determined from the sum of capacitances through the path 
is the lowest. Similar results were previously reported in CPP measurements 
for nonmagnetic granular films [93]. It is considered that at higher voltages 
the contribution of various current paths including those with higher charging 
energy appears to lead to the irregular period of Coulomb staircase. 

Fig. 5-8 (b) shows the bias voltage dependence of the TMR derived from 
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the two I — Vb curves shown in Fig. 5-8 (a). We observed the oscillation of the 
TMR as a function of bias voltage. The peak of TMR repeats with the period 
of the Coulomb staircase. The modification of the I — Vb curve by the applied 
field seems to bring about the enhancement of the TMR with the steps of the 
Coulomb staircase, resulting in the oscillation of the TMR. The largest TMR 
(about 20 %) was seen at Vb = 15 mV; however, TMR converged to almost zero 
as the bias voltage was further increased. This is probably because the barrier 
quality of the Al-0 bottleneck layer is still poor, leading to the rapid decrease in 
TMR with bias voltage. It is also noted that a sign change was observed in the 
TMR. One possible origin we may consider is the effect of spin accumulation. 
In order to clarify TMR behavior, we performed a further study considering the 
spin accumulation effect in a nanoparticle. 



24 



6. Spin accumulation effect in nanoparticles 

We observed TMR oscillation associated with the Coulomb staircase in CPP 
geometry samples fabricated by FIB. In this chapter, we report experimental 
evidence for the spin accumulation effect in Co nanoparticles leading to the 
oscillation of TMR with alternate sign changes. Furthermore, the spin relax- 
ation time in Co nanoparticles was evaluated by comparing the experimental 
results to numerical simulations. We discovered that the spin relaxation time 
in the nanoparticles is unprecedentedly enhanced up to more than the order of 
hundreds of nano seconds, compared to that evaluated from the spin diffusion 
length of ferromagnetic layers in previous CPP-GMR studies, i.e., the order of 
tens of picoseconds. 

6.1 Experimental procedures and results 

In order to study this subject, the sample design and the fabrication pro- 
cess were modified. We fabricated pillared structures consisting of Al electrode 
/ Al-0 / Co-Al-0 granular film / Co electrode layers (Fig. 6-1) to measure 
their magnctotransport properties in CPP geometry. Samples were prepared on 
thermally oxidized silicon substrates as follows: the bottom Al electrode was 
first deposited using ion beam sputtering (IBS), and the surface was plasma 
oxidized to form a conductance bottleneck layer for the observation of Coulomb 
staircase. A 15 nm thick Co-Al-0 granular film was then deposited by reac- 
tive rf sputtering. The top Co electrode and Pt layer were deposited by IBS. 
This layered structure was then microfabricatcd to reduce the contact area us- 
ing electron-beam (EB) lithography and an Ar ion etching process. An Al / 
Al-0 / Co-Al-0 / Co / Pt pillar 0.4 fj,m x 0.4 ^im in area was prepared by 
the following process: EB patterning on positive resists, depositing of etching 
mask (Co), lift-off and Ar ion etching. To prevent short-circuiting between the 
electrodes, the sidewall of the pillar was then plasma-oxidized and covered with 
a thick Al-0 also deposited by IBS. Finally a Pt layer was deposited to contact 
the top electrode. The sample structure is schematically shown in Fig. 6-1. 

Fig. 6-2 (a) shows the I — ^characteristics at H— (black line) and H= 
10 kOe (gray line) in the positive bias voltage region. The / — ^curves reveal a 
definite Coulomb staircase, which was also observed in the negative bias voltage 
region (not shown here). The first step of the Coulomb staircase appears at 
about 15 mV, and the subsequent steps appear at 50 mV, 85 mV and 120 mV, 
respectively. We think the current preferentially flowed through a single or a 
few restricted local paths where only one particle was located between electrodes 
for the following reason: When the nominal thickness of the Co-Al-0 layer was 
smaller than 15 nm, the samples show neither Coulomb blockade nor Coulomb 
staircase. In other words, 15 nm is the minimum thickness to observe the single- 
electron phenomena. This suggests that the Co-Al-0 layer has considerable 
thickness fluctuation, and therefore there are some thin parts where only one 
particle exists in the direction normal to the film plane. (If the thickness of 
the Co-Al-0 layer is completely uniform, three or four particles should exist in 
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the direction normal to the film plane everywhere, and the critical thickness for 
SET should be much smaller than 15 nm because the particle size is a few nm.) 
Furthermore, the Al-0 bottleneck layer also has thickness fluctuation, playing 
a role as a path restriction filter. The shape of the staircase measured at H= 
is different from that of at H— 10 kOc. Comparing the shapes of these two 
curves, it appears that the steps of the staircase at H= are steeper than those 
at H= 10 kOe, leading to a periodic intersection of the two curves around the 
step points of the staircase. 

Fig. 6-2 (b) shows the ^dependence of TMR derived from the I — Vcurves 
shown in Fig. 6-2 (a). TMR oscillates with the same period as the staircase, 
ranges from -10 % to +15 % and shows alternate sign changes: negative peaks, 
i.e., inverse TMR, appear around the step points of the staircase, whereas pos- 
itive bumps appear between these negative peaks. In a conventional ferromag- 
netic tunnel junction of macroscopic size, TMR should be limited within that 
expected from the spin-polarization factors of ferromagnetic components; and 
it decreases monotonically with bias voltage. The expected value of the TMR 
in this sample should be less than 2% because the resistance of the bottleneck 
which is located at the spin- independent junction between the Al electrode and 
a Co particle is 10 times larger than that of the other spin- dependent junction. 
Nevertheless, TMR ranging from -10 % to +15 % was observed, indicating a 
significant enhancement and anomalous oscillation with sign change of TMR. In 
order to confirm the sign of TMR, the magnetic field dependence of the electri- 
cal resistance (MR curves) was measured. Fig. 6-2 (c) shows MR curves at V= 
0.05 V and 0.12 V at which TMR shows negative and positive values, respec- 
tively. Both MR curves show a hysteretic behavior with two peaks according 
to the coercive force of the Co nanoparticles. The MR curve at V= 0.05 V 
shows an increase of electrical resistance (R) with increasing applied magnetic 
field; i.e., inverse TMR is observed, whereas the MR curve at V= 0.12 V shows 
a decrease of i?with the magnetic field. It is noted that the shift of MR curves 
toward the positive direction of the magnetic field is seen. A possible origin for 
the shift of the MR curves is the exchange bias between Co and CoO formed by 
the surface oxidization of the Co nanoparticles. 

6.2 Numerical analysis 

We performed numerical calculations on the basis of the orthodox theory 
with an arbitrary spin relaxation time in ferromagnetic tunnel junctions and 
compared it with the ^dependence of TMR in the present experiment. Since in 
the conducting region the Co- Al-0 film contains only one nanoparticle between 
electrodes as mentioned before, we employed a PM / I / FM particle / I / FM 
double-tunnel-junction model (see Fig. 6-3 (a)) to analyze the experimental 
data. We considered the collincar alignment of the magnetizations. Although 
Fig. 6-2 (c) shows that the angle between the magnetizations of the Co particle 
and the Co electrode is large at zero magnetic field, it is hard to conclude that 
the magnetizations were aligned to be anti-parallel since the easy magnetiza- 
tion directions of nanoparticles were distributed randomly and the direction of 
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the applied magnetic field was fixed to a certain direction. If the magnetiza- 
tions are non-collinearly aligned the tunneling resistance for each spin depends 
on the relative angle of magnetizations [11, 94] and we have to calculate the 
nonequilibrium spin distribution matrix [95, 96]. However, it is important to 
note that the calculations for collinear configurations give the lower limit of the 
spin relaxation time. Hereafter, we consider only the collinear configurations. 

Without losing generality we can take the spin quantization axis so that the 
minority and majority spins of the nanoparticle are represented by J. and j, 
respectively. The spin accumulation occurs in the ferromagnetic island where 
there is a significant difference in the densities of states for minority (-Dj.) and 
majority (D-p) spins at the Fermi level due to the exchange splitting in the 
c?-band. According to the band calculation of Co, the net spin polarization 
Pdos = [D] — D^)/(D^ +Di) is taken to be -0.73. In addition to Pdos > another 
spin polarization must be taken into account, i.e., that for tunneling electrons 
which is dominated by s — ^hybridization [97]. It is described as Ptun, hereafter. 
PtunOi Co is assumed to be 0.35 [98]. We suppose that the energy relaxation 
time is so short that the Fermi distribution function can be used to represent the 
electron distribution in the Co nanoparticle. The model shown schematically in 
Fig. 6-3 (a) assumes that the magnetization vector of the nanoparticle is in a 
fixed orientation, whereas that of the right electrode is reversed by applying a 
magnetic field, leading to parallel (P) and anti-parallel (AP) alignments. The 
TMR is defined as TMR = 1 - R P /R AP , where R P {AP) is the resistance of the 
whole system in the P (AP) alignment. The tunnel resistance and capacitance 
of each junction have been adjusted to reproduce the magnitude of the tunneling 
current and the period of the Coulomb staircase. The capacitances are taken 
to be Ci = 4.44 aF and Ci = 3.00 aF, where subscripts 1 and 2 correspond to 
the left and right junctions, respectively. The tunnel resistances for the parallel 
alignment were determined in the following way: first we decided that the tunnel 
resistance of the majority spin for parallel alignment at the left junction (-Rf •)■) 
was 25.0 GO, and the conductance bottleneck was placed at the left junction 
(PM / I / FM particle), assuming that i?fj was 10 times larger than R^- The 
other resistances are given by 

Rl i =Rl i {l + PAl){l + Ppartcle)/0--PAl)0--Ppartcle) = 51.9GO, 

i?£ T =0.1 -i?f T = 250MO, 

R^= 0.1 ■R^(l+PparUcle)(l + Pco)/ (I -PparUcle)(l-Pco) = 1.08Gfi, 

(6.1) 

where Pco, Pparticle an d Pai are the Ptun of the right electrode, the parti- 
cle and the left electrode, and are assumed to be 0.35, 0.35 and 0, respectively. 
These considerably large resistances (GO) are reasonable for such a small tunnel 
junction consisting of nanoparticles. For the antiparallel alignment, tunnel resis- 
tances at the left junction were Rf 1 ^ — i?f ^ = 25.0 GO and R^ = i?f ^ = 51.9 
GO. At the right junction, D Partic i ed = Dp arWc(e>T (1 - Pp ar tcile) I (1 + Pparticle), 
D Co ,i = Dco,t (1 - Pco) I (1 + Pco) and R^ cx (L>p Qrtlc z e , T L> C 'o,T)~ 1 produced 
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the following resistances: 



R$A =R^(l + Pco)/(l-Pco) = 519MQ, 

Rtl = R 2,f( 1 + P Particle)/(l-PparUcle) = 519MO. (6.2) 

When an external voltage is applied, electrons flow from the left to the right 
electrode. The current flowing in the junction is calculated using the master 
equation technique [21, 22] . By applying the charge conservation and Kirchhoff 's 
law and by using the Fermi golden rule for the calculation of the tunneling 
Hamiltonian, we could write the tunneling rate for the i-th junction as follows: 

l ' A ' e 2 iV exp {Ef(n) + AE F /k B T) 1 { ' 

where i stands for the considered junction (i= 1, 2), a stands for the direction of 
the spin, and the exponent ± indicates whether an electron is added or removed 
from the particle. Finally, the energy is written: 

Ef{n) = (1 + 2n)^ ± ^eV (where j ± i) (6.4) 

where C = C\ + C*2, T is the temperature, and AE F the shift of the Fermi 
energy of the island due to spin accumulation for the cr-spin population. In 
the stationary regime, the transition rates for incoming and outgoing electrons 
cancel each other thes and the probability p(n) of finding n additional electrons 
in the central particle may then be calculated. Finally, the current is expressed 
as follows: 

oo 

J=-e£ E (It»-ri;») (6.5) 

a n— — oo 

The spin accumulation in a Co nanoparticle is represented by the spin-dependent 
shift of the chemical potential for thegergpin state AE F , which is determined 
by charge-neutrality and spin-conservation conditions. The former is given by 
D^AEp = —D^AEp, where D a is the density of states at the Fermi level for 
spin a. The latter is expressed as follows: 

~ = ^AE° F (6.6) 
e t sf 

where Ii^ a is the current at the i th junction (i= 1, 2) for spin a, is the volume 
of the particle, tsf is the mean spin relaxation time in a nanoparticle defined as 
Tg F = (t^ 1 +Tj^ 1 )/2, and r a is the spin relaxation time of electrons with spin a. 
For a ferromagnet, tj and t± are not the same but satisfy the detailed balancing 
equation, D^/t^ = DJt^. The tunnel current is obtained by solving the master 
equation, charge-neutrality and spin-conservation conditions, self-consistently 
[14, 15, 27, 36]. 

We first show the results of numerical calculation assuming an infinite spin 
relaxation time in Fig. 6-3 (b), (c) and (d). The significant effect of the spin 



28 



accumulation is found in these results. Fig. 6-3 (b) shows the Independence 
of the chemical potential shifts AE F for the majority and minority spins in 
the P and AP alignment. It shows a sawlike oscillation with a period close 
to that of the Coulomb staircase shown in Fig. 6-3 (c). The important point 
is that the difference in the energy shift AE F for each alignment and spin 
leads to a splitting of the discrete charging levels Ec' Pi " AP \n) for spin a and 
P(AP) alignments, where ndenotes the numbers of the charge in the particle. 

j^a,P(AP) ^ . g CX p ressec [ ag 

EV p{A P){n) = (n+i/iy + AE « F P(AP) {n (67) 

Cl 

The voltage corresponding to the step points of the Coulomb staircase is then 
accordingly split. In Fig. 6-3 (c), one can see that each I— V curve has a unique 
shape because AE F affects the tunnel conductance, and furthermore intersects 
around the step points of the staircase due to the split of EZ' P ^ AP \ri). This in 
turn leads to the periodical sign change of the ^dependence of TMR (TMR- 
Vcurve) with a period close to that of the Coulomb staircase as shown in Fig. 
6-3 (d). 

We also performed a calculation assuming an infinite spin relaxation time in 
the case of the converse arrangement of the conductance bottleneck; i.e., i? 2 j 
was fixed to be 25.0 Gf2 and was assumed to be 10 times larger than R p The 
other resistances and capacitances are described in Fig. 6-4 (a). We show the 
results in Fig. 6-4 (b), (c) and (d). Though the V dependence of the chemi- 
cal potential also appears as sawlike oscillations with a period of the Coulomb 
staircase, none of the I—V curve intersect; thus no sign changes occur for TMR 
in the entire V range. These results indicate that the alternate sign change of 
TMR due to spin accumulation in the Co nanoparticle is reproduced when the 
conductance bottleneck is put at the junction of the PM / I / FM particle. 

6.3 Evaluation of spin relaxation time 

Fig. 6-5 shows the TMR-Vcurve for various values of the spin relaxation 
time tsf- Although the period of the oscillation does not depend on tsf, the 
shape of the TMR-Vcurve is quite different between the fast and slow spin 
relaxation regimes. The critical value of the spin relaxation time which divides 
the fast- and slow-spin relaxation regimes is about 100 ns, which corresponds to 
the interval between two successive tunneling events (e/I~ 100 ns). In the fast- 
spin relaxation regime (tsf < 100 ns), the TMR-Vcurve shows sharp positive 
peaks at the step points of the Coulomb staircase. On the contrary, in the slow 
spin relaxation regime {tsf > 100 ns), sharp negative peaks appear at the 
same step points of the staircase, and broad positive bumps appear between 
two successive negative peaks. In Fig. 6-6 (a), (b), (c) and (d), we show cross- 
sections of Fig. 6-5 at tsf = 1 ns, 10 ns, 150 ns and infinity, respectively. At 
tsf = 1 ns, in the fast-spin relaxation regime, we have sharp positive peaks at 
V= 18,53,88, and 123 mV, where the tunneling current shows steps. The TMR- 
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V curve maintains the same shape as far as the spin relaxation time tsf < 1 ns. 
However, at t$f = 10 ns, one can see that dips which take negative TMR values 
appear at the step points of the Coulomb staircase. At t$f — 150 ns, the TMR 
shows alternate sign changes with moderate amplitude reaching a maximum 
at infinite tsf- The negative peaks associated with the steps of the Coulomb 
staircase provide clear evidence of spin accumulation in Co nanoparticles [15, 
27]. Finally as shown in Fig. 6-7, the experimentally observed TMR curve is 
well reproduced by choosing the spin relaxation time to be tsf — 150 ns. 

The present TMR experiments have shown that the tunnel current through 
nanoparticles is drastically affected by spin accumulation, and have revealed 
that the TMR exhibits a sign change at the step points of the Coulomb stair- 
case. Although the spin relaxation mechanism of the magnetic nanoparticle is 
not well understood at present, there are some candidates for the cause of such a 
dramatic increase of the spin relaxation time. One is the quantization of energy 
levels, which is the most characteristic feature of nanoparticles with diameters of 
1^5nm [99-102]. The discrete energy level spacing in Co nanoparticles having a 
mean diameter of ^2.5 nm is estimated to be about 2 meV, which is larger than 
the thermal fluctuation at 4.2 K (<~0.4 meV). The discreteness of the energy lev- 
els due to the zero dimensionality of nanoparticles leads to strong suppression 
of spin-flip scattering by the spin-orbit interaction in the nanoparticles [103, 
104]. Mitrikas et al. have studied the spin-lattice relaxation of paramagnetic 
nanoparticles embedded in amorphous (Si02) and a crystalline (Ti02) matrices 
[105]. They showed that the spin lattice relaxation (SLR) is blocked due to the 
amorphous phase of the matrix. Although their nanoparticle is not ferromag- 
netic, the blocking of the SLR is another candidate for the strong enhancement 
of the spin relaxation time in our sample because the insulating Al-0 matrix 
is in an amorphous phase. Furthermore, the spin-flip mechanisms due to the 
interaction with magnons in a nanoparticle are also suppressed because of the 
size quantization of magnon excitations [106]. 

Recently, further experimental studies on spin accumulation in paramagnetic 
Au nanoparticles have been performed [68, 72]. Observation of non-zero TMR 
is attributed to the spin accumulation in Au nanoparticles in an FM / I / Au 
particles / I / FM structure, which should not show TMR in the framework 
of Jullicre model. Although these studies did not mention the spin relaxation 
time in Au nanoparticles, the nano-ampere order of the observed current sug- 
gested a nano-second order for the spin relaxation time (e/J). It is suggested 
to be enhanced even in the Au nanoparticles. All of these particular features of 
nanoparticles can lead to an extremely long spin relaxation time and a striking 
spin accumulation effect on TMR. To determine the dominant mechanism of 
spin relaxation in nanoparticles, further studies investigating the dependence 
of spin relaxation time on variables such as the temperature, the material of 
the insulating matrix and particle size are needed. We note that the results 
shown here are promising for the potential application of nanoparticles as basic 
elements of spin-electronic devices. 
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7. Summary 

We have reviewed studies on spin-dependent transport in systems contain- 
ing ferromagnetic nanoparticles. In a tunnel junction with a ferromagnetic 
nano-island and electrode, spin-dependent single-electron tunneling (SD-SET) 
gives rise to remarkable tunnel magnetoresistance (TMR) phenomena. We stud- 
ied magnetotransport properties in both sequential tunneling and cotunncling 
regimes of SET and found the enhancement and oscillation of TMR. 

The self-assembled ferromagnetic nanoparticles we used in this study con- 
sisted of a Co-Al-0 granular film, with cobalt nanoparticles embedded in the 
Al-0 insulating matrix. In Ch. 3 we showed the preparation methods, and basic 
transport and magnetic properties of the films. A C036AI22O42 film prepared 
by the reactive sputtering method produced a TMR ratio reaching 10 % and 
superparamagnetic behavior at room temperature. The fitting of the magne- 
tization curve to that of the calculation from the Langevin function revealed 
the distribution of particle sizes ranging from 2 nm to 5 nm. The TMR ratio 
exhibited an anomalous increase at low temperatures but no significant indi- 
cation of change with bias voltage. In Ch. 4, we showed that the anomalous 
increase of the MR indicated evidence for higher-order tunneling (cotunncling) 
between large granules through intervening small granules. We emphasize that 
the existence of higher-order tunneling is a natural consequence of the granular 
structure, since broad distribution of granule sizes is an intrinsic property of 
granular systems. 

In Ch. 5, we concentrated on SD-SET properties in a sequential tunneling 
regime. In order to study it, the tunnel paths between electrodes had to be re- 
stricted because the film contains a large number of particles of different sizes. 
We fabricated two types of device structures with Co- Al-0 film using focused 
ion-beam milling or electron-beam lithography techniques. One is a granular 
nanobridge structure: point-shaped electrodes separated by a very narrow lat- 
eral gap filled with Co- Al-0 granular film. The other is a current-perpendicular- 
to-plane (CPP) geometry structure: a thin Co-Al-0 granular film sandwiched 
by ferromagnetic electrodes with the current flowing in the direction perpendic- 
ular to the film plane through a few Co particles. We measured the current-bias 
voltage curves in these samples, and found the enhancement and oscillation of 
TMR due to spin-dependent SET in a sequential tunneling regime. We made 
a theoretical calculation employing the orthodox theory of SET and explained 
the enhanced TMR just above the Coulomb threshold voltage. In Ch. 6, we 
reported experimental evidence of the spin accumulation effect in Co nanoparti- 
cles leading to the oscillation of TMR with alternate sign changes. Furthermore, 
the spin relaxation time in Co nanoparticles was also evaluated by comparing 
the experimental results to numerical simulations. We discovered that the spin 
relaxation time in the nanoparticles is unprecedentedly enhanced up to more 
than the order of hundreds of nanoseconds, compared to that evaluated from 
the spin diffusion length of ferromagnetic layers in previous CPP-GMR studies, 
i.e., the order of tens of picoseconds. 

Although the spin relaxation mechanism of the magnetic nanoparticle is 
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not well understood at present, there are some candidates for cause of such a 
dramatic increase of the spin relaxation time: suppression of spin-flip scatter- 
ing by spin-orbit interaction due to the discrete energy level spacing in the Co 
nanoparticle and/or blockade of spin-lattice relaxation due to the amorphous 
Al-0 matrix. Further studies such as measuring the temperature dependence of 
spin relaxation time will clarify the mechanism of spin relaxation in nanoparti- 
cles. 
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Figure captions 



Fig. 2-1 Schematic diagram of a circuit including a single-electron tunneling 
(SET) double junction. 

Fig. 2-2 Current (I) - bias voltage (V") characteristics of a SET junction (a) 
with no asymmetry [of tunnel] resistances, and (b) with strong asymmetry of 
tunnel resistances. 

Fig. 2-3 Schematic diagram of a double tunnel junction to explain chemical 
potential shift due to spin accumulation in the ferromagnetic nanoparticle. 

Fig. 3-1 Schematic illustration of a sputtering method for preparing insu- 
lating granular films: (a) reactive sputtering, (b) sputtering with a composition 
target and (c) tandem deposition with plural targets. 

Fig. 3-2 (a) Plan view and (b) cross-sectional view of transmission electron 
microscopy (TEM) micrographs for a C046AI19O35 film. 

Fig. 3-3 High-resolution electron microscopy (HREM) micrograph for Co 36 Al 2 2 
granular film. 

Fig. 3-4 Temperature (T) dependence of electrical resistivity (p) for Co-Al-0 
films of different compositions: C054AI21O25, C052AI20O28, C046AI19O35, and 
C036AI22O42. 

Fig. 3-5 (a) Magnetic-field dependence of magnetization (M — H) and (b) 
temperature dependence of magnetization (M — T) for the C036AI22O42 film. 
In (b), M — T curve of field cool (FC) was measured at H — 20 Oe. 

Fig. 3-6 Magnetic-field dependence of electrical resistivity (MR curve) mea- 
sured at (a) room temperature and (b) 4.2 K. Magnetoresistance (MR) ratio 
shown in the left axis was estimated from the definition: MR = Apvp max . 

Fig. 3-7 The result of particle size distribution of a Co46Alig0 3 5 film esti- 
mated from superparamagnetic behavior, (a) Result of fitting calculated mag- 
netization curves to experimental results at T — 200 and 293 K. (b) The size 
distribution which gives the best fit in (a). 

Fig. 4-1 Temperature (T) dependence of the tunnel magnetoresistance (TMR) 
ratio for Co-Al-0 films of various compositions: Co 54 Al2i0 2 5, C052AI20O28, 
Co4 6 Al 19 03 5 , and Co 36 Al 2 2 42. Solid curves (a, b, c, and d) represent the the- 
oretical magnetoresistance ratio given by Eq. (4.2) with spin polarization P; a: 
P = 0.306; b: P = 0.290; c: P = 0.275; d: P = 0.250. 
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Fig. 4-2 Schematic illustrations of current-perpendicular-to-plane (CPP) ge- 
ometry sample where a 1 ^m thick C036AI22O42 granular film was sandwiched 
between upper and lower Au-Cr electrodes. 

Fig. 4-3 Bias voltage (Vb) dependence of p at H = (a) and V& dependence 
of TMR (b) for a CPP geometry sample with C036AI22O42 measured at 4.2 K. 
Closed circles represent the experimental results and solid curves represent the 
theoretical ones using Eq. (4.6) withT = 4.2 K and N g = 140. 

Fig. 4-4 (a) Schematic illustration of granular structure and a higher-order 
tunneling process where a charge carrier is transferred from the charged large 
granule (left), via the two small ones, to the neutral large one (right), (b) Model 
structure used for the calculation of conductivity 

Fig. 4-5 Temperature (T) dependence of electrical resistivity (p) for Co-Al-0 
films of different compositions: C054AI21O25, C052AI20O28, C046AI19O35, and 
C036AI22O42. Solid lines represent the theoretical results given by Eq. (4.2) 
with charging energy E c ; Co 36 Al 2 2 42: E c / k B = 110 K; Co46Ali 9 35 :£; c / 
k B = 25 K; Co 5 2Al 20 O2 8 : E c / k B = 18 K; Co 54 Al 21 2 5: E c / k B = 9 K. 

Fig. 5-1 (a) Scanning tunneling microscopy (STM) topographic image and 
(b) current - bias voltage characteristics based on scanning tunneling spec- 
troscopy (STS) measurements at room temperature for a C036AI22O42 film. 

Fig. 5-2 (a) Schematic view of an insulating granular nanobridge. (b) Scan- 
ning ion microscopy image of the NbZrSi electrodes separated by a nanometer- 
sized lateral gap. The length (I) is 30 nm and the width (w) is 60 nm. 

Fig. 5-3 (a) Current - bias voltage (I — Vb) characteristics and (b) H de- 
pendence of TMR measured at 4.2 K for a sample of w = 60 nm, I = 30 nm 
and t — 7.5 nm. In (a), the solid and dashed curves represent I — V> curves in 
a magnified current range at H = and 10 kOe, respectively. I — V> curves 
throughout the measured current range are shown in the inset. 

Fig. 5-4 (a) Current - bias voltage (7 — Vb) characteristics and (b) Vj, de- 
pendence of TMR measured at 4.2 K for a sample of w = 700 nm, I = 40 nm 
and t = 15 nm. Vicwgraphs are shown in the same form as in Fig. 5-3. 

Fig. 5-5 V p vs. V t h for the various samples. V p is the voltage where the 
TMR shows the maximum. 

Fig. 5-6 (a) Schematic illustration and equivalent circuit of the calculation 
model. The conducting path in the contact is modeled by a parallel circuit of 20 
triple-tunnel-junctions, assuming the tunnel resistance between an electrode and 
a particle (R ep ), tunnel resistances between particles (R pp ) and junction capac- 
itances (C ep and C pp ). These parameters are assumed to be distributed around 
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the mean values, i.e., R ep — (l±5)R ep , R pp — (1±5)R PP , C ep — (l±5)C ep , and 
Cpp = (1± S)C PP , where 5 is the deviation from the typical value, (b) The nu- 
merical calculation results (solid curve) and experimental results (dashed curve). 

Fig. 5-7 Schematic illustrations of (a) the CPP geometry sample, and (b) 
the process of making a contact area by FIB etching. 

Fig. 5-8 (a) Current - bias voltage (I — Vb) curves and (b) bias-voltage 
dependence of TMR measured at 4.2 K for a CPP geometry sample with a 0.5 
x 0.5 /im 2 contact area. In (a), the solid and dashed curves represent I — Vb 
curves at H = and 10 kOe, respectively. 

Fig. 6-1 Schematic illustration of a sample with pillar structure, prepared 
for the current-perpendicular-to-plane (CPP) measurement. It consists of Al 
bottom electrode / Al-0 / Co-Al-0 granular film / Co top electrode / Pt mi- 
crofabricated by electron-beam lithography and Ar ion milling process to reduce 
the contact area down to 0.4 x 0.4 ^m. 

Fig. 6-2 Magnetotransport properties measured at 4.2 K in an Al / Al-0 / 
Co-Al-0 / Co sample of sub-micron sized area, (a) Current (I) - bias voltage 
(V) curves and (b) V dependence of the tunnel magnetoresistance (TMR). In 
(a), the gray and black solid curves represent a clear Coulomb staircase at zero 
magnetic field (H = 0) and applied field (H = 10 kOe), respectively, (c) The 
magnetic field dependence of electrical resistance (MR curves) at V = 0.05 V 
and 0.12 V at which TMR shows negative and positive values, respectively. 

Fig. 6-3 Numerical calculation results considering the spin accumulation 
effect on TMR. The model assumes that a spin relaxation time in the ferromag- 
netic island is infinite and that a conductance bottleneck exists at the junction 
between the ferromagnetic island and the paramagnet electrode, (a) Schematic 
illustration of the model: paramagnet / ferromagnet / ferromagnet double tun- 
nel junction, (b) V dependence of the chemical potential splitting AE F , which 
is caused by the spin accumulation, (c) Current (I) - bias voltage (V) curves in 
the anti-parallel (AP, dashed line) and parallel (P, solid line) alignments, (d) V 
dependence of the tunnel magnetoresistance (TMR). 

Fig. 6-4 Numerical calculation results in the case of converse arrangement 
of the conductance bottleneck, i.e., R^P is fixed to be 25.0 Gf2 and to be 10 
times larger than i?f (a) Schematic illustration of the model: paramagnet 
/ ferromagnet / ferromagnet double tunnel junction, (b) V dependence of the 
chemical potential spliting AE F . (c) Current (I) - bias voltage (V) curves in 
the anti-parallel (AP, dashed line) and parallel (P, solid line) alignments, (d) V 
dependence of TMR. 

Fig. 6-5 Surface plot of the TMR, which is calculated by the orthodox the- 
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ory, as a function of V and tsf- 

Fig. 6-6 V dependence of TMR for various values of the spin relaxation time 
tsf- (a) 1 ns, (b) 10 ns, (c) 150 ns and (d) infinity. 

Fig. 6-7 V dependence of TMR of [the] calculation result (tsf — 150 ns) 
and experimental TMR. The TMR curve is well reproduced by choosing the 
spin relaxation time to be tsf = 150 ns. 
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